X-ray afterglow light curves : toward standard candle ? 
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ABSTRACT 

We investigate the clustering of afterglow light curves observed at X-ray and optical wave- 
lengths. We have constructed a sample of 61 bursts with known distance and X-ray afterglow. 
This sample includes bursts observed by BeppoSAX, XMM-Newton, Chandra, and SWIFT. We 
correct the light curves for cosmological effects and compare the observed X-ray fluxes one day 
after the burst. We check for correlations between the observed flux and the burst spectral and 
temporal properties. We confirm the previous result of Boer and Gendre (2000) that X-ray after- 
glow light curves cluster in luminosity, even when we consider the last SWIFT data. We observe 
this clustering only for the afterglow light curves; the inclusion of prompt-related data broaden 
the distribution. A similar clustering is observed for the optical light curves; GRB sources can be 
divided in three classes, namely optical and X-ray bright afterglows, optical and X-ray dim ones, 
and optically bright -X-ray dim ones. We argue that this clustering is related to the fireball total 
energy, the external medium density, the fraction of fireball energy going in relativistic electrons 
and magnetic fields. These parameters can be either fixed to a standard value, or correlated. We 
finally propose a method for the estimation of the GRB source redshift based on the observed 
X-ray flux one day after the burst and optical properties. Using this method, we compute a 
redshift of 1.4 ± 0.2 for GRB 980519 and of 1.9 ± 0.3 for GRB 040827. We tested this method 
on three recently detected SWIFT GRBs with known redshift, and found it in good agreement 
with the reported distance from optical spectroscopy. 

Subject headings: gamma-ray:bursts - X-ray:general 



1. Introduction 

Long Gamma-Ray Bursts (GRBs) are linked 
wit h the death of massive stars (for a review, 
see Meszaros 2006). Their association with su- 



perno vae (e.g. Hiorth et al. 2003t IStanek et al. 
20031 ) and the fact that these events are at cosmo- 
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logical distance (|Metzger et al. 1119971) make them 
interesting for studies of cosmology in the red- 
shift range 1-15. However, while the GRB fea- 
tures a prompt emission, usually seen in gamma- 
ray only, and an afterglow seen at all wavelengths, 
only the former emission has been considered for 
cosmological studies yet. One of the first at- 
tempt to do so, based on the — E^^n corre- 



lat ion found by lAmati et al. I (j2002i ). was done 



by iGhirlanda et al. I (|2004l ) , who constrained the 
fim — f^A parameters. While not very constraining, 
their findings were compatible with previous tests 
made with supernova samples. The use of GRB 
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for cosmological studies, specially at high redshift 
where their detection at optical/IR wavelength is 
difficult, needs to build a robust indicator of their 
distance, whenever possible based on their intrin- 
sic properties. 

The afterglow emission is less studied from a 
cosmological point of view because of its diver- 
sity. However, there were hints of standardiza- 
tion of the X-ray afterglows. The first attempt 
was done by iBoer fc Gendre ( 20001 ) , who found 
evidences for clustering in the X-ray light curves 
of BeppoSAX after glows. This clusterin g was 



confirmed later by Gendre fc Boer ( 2005^ who 



extended the sample to the XMM-Newton and 
Chandra data. In the following we will refer to 
these a rticles respectively as pape r I and II. In ad- 
dition, iKouve liotou et al. I (j2004l ) showed that su- 
pernova and GRB light curves had similar behav- 
iors and were converging with time toward a sim- 
ilar luminosity. In paper I we tried also to check 
if optical light curves were also clustered: this at- 
tempt failed due to the poor knowledge of the in- 
trinsic absorption in the burst host gala xy at that 



time. This study was completed by Na rdini et al. 



(|2006l ) and iLiang fc Zhimg] (|2006f) who found in 
dependently that optical afterglows were also clus- 
tered in luminosity. 

The high detection rate and throughput of 
SWIFT, which provide a large sample of X-ray 
observations and rapid, accurate localization, en- 
abling a redshift estimation by optical/IR tele- 
scopes, allowed us to increase dramatically our 
sample. With this larger sample, we have tried 
to derive a method for estimati ng the burst red- 



shift from the X-ray light curve ([Gendre fc Boer 
[200a . However, before to apply this method it is 
necessary to understand the nature of these two 
groups, and to determine a consistent way to de- 
rive to which group belongs each burst. This is the 
purpose of this article. In Section[2]we present our 
sample and the data analysis we performed. We 
discuss the X-ray clustering in Section [S] We com- 
pare the X-ray and optical clusterings in Section 
m In section [5] we present our method of GRB 
source redshift estimation from the X-ray after- 
glow light curve. We finally test this method on 
several GRBs detected recently by SWIFT, and 
we propose an estimation of the redshift for two 
previously detected GRBs of unknown distance. 



2. X-ray afterglov^r sample and analysis 

Our sample of GRBs with known redshift and 
X-ray afterglow observations is listed in Table 
[TJ It includes all afterglows observed by Bep- 
poSAX, XMM-Newto n and Chandra, whi c h data 
were 



retrieved fro m 
and iGendre et al. 



.De Pasquale et al. 
(2006|). As for SWIFT 



I 200d ) 

ob- 



servations, we browsed the SWIFT archive web 
page and selected all bursts detected prior to the 
1st of August 2006 with a measured redshift and 
a igo larger than 2.0 seconds (in order to exclude 
short bursts). From this sample of SWIFT bursts, 
we further excluded GRB 060123 due to data pro- 
cessing errors. We removed the flaring parts of 
the light curve when applicable (this totally re- 
moved from the sample the data of GRB 050904, 
which is thus not listed in Table [T]). We corrected 
the flux light curves for distance effects as in Pa- 
per I and II : we apply a k-correction using the 
measured spectral index of each afterglow, assum- 
ing a flat universe with = 0.3, and correcting 
for the time dilation effect by computing the lu- 
minosity at a time t using the observed flux at 
time t X (1 + z). This is very important for the 
correction accuracy, as discussed in Sec. U) We 
restricted the light curves to the 2.0—10.0 keV 
X-ray band, where the absorption is negligible. 
This allowed us to neglect any other corrections 
for absorption by the ISM. We do not take into 
account any beaming effect due to a possible jet 
structure. 

As we done previously, instead of using a lu- 
minosity light curve, we express all light curves 
in flux units at a given distance (like the optical 
absolute magnitude). This allows to reduce the 
uncertainties on the correction. For consistency 
with paper I and II, we fix the redshift to z = 1. 
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Table 1 

The burst sample used in this paper. We indicate for each burst the satellite that 
observed the x-ray afterglow, its redshift, spectral and decay index, its group 



assignation (see text ) IN X-RAYS AND IN OPTICAL (EXTRACTED FR OM INARDINI ET AL. 112006 



Liang fc Zhang Ibood ). and the values of Eiso and Ep listed in IAmati^ ( 2006[) . For the 



BURSTS DETECTED BY SWIFT WE ALSO INDICATE THE Ta VALUES LISTED IN I WILLING ALE ET AL. 



(j2007l ) CORRESPONDING TO THE AFTERGLOW STAR T . THE SPECTRAL AN D DECAY INDEX IN F ORMA TION 



ARE EXTRACTED FROM 



De Pasquale ET AL. I (|2006|. BeppoSAX). IGeNDRE ET AL. I (|2006 



XMM-Newton. Chandra! , and IWillingale et al~I (j2007l . SWIFT). 



Burst name 


X-ra,v 


redshift 


Decay 


Spectral 


X-ra,v 


Optical 


E 




los(T„ ) 




satellite 




index 


index 




group 


group 


(10^^ erg) 


(keV) 


(s) 


GRB 970228 


BeppoSAX 


0.695 


1.3 ± 0.2 


1.0 ±0 


3 


II 


I 


1.84 ±0.14 


195 ± 64 





GRB 970508 


BeppoSAX 


0.835 


0.8 ±0.2 


1.4 ±0 


3 


II 


II ? 


0.71 ± 0.15 


145 ± 43 





GRB 971214 


BeppoSAX 


3.42 


1.0 ± 0.3 


1.1 ± 
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685 ± 133 
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BeppoSAX 


1.60 


1 45-1-0 06 


0.99 ± 


05 


I 


I 




1 724 + 466 




GRB 990510 


BeppoSAX 


1.619 


1.4 ± 0.1 


1.17±0 


09 


I 


I 


20 ± 3 


423 ± 42 




GRB 991216 


Chandra 


1.02 


(0.9) 


0.7 ±0 


2 


I 


I 


78 ±8 


648 ± 134 





GRB 000210 


BeppoSAX 


0.846 


1.38 ±0.03 


0.9 ±0 


2 


II 





17.3 ± 1.9 


753 ± 26 





GRB 000214 


BeppoSAX 


0.42 


0.7 ± 0.4 


1.0 ±0 


3 


II 





0.93 ± 0.03 


> 117 





GRB 000926 


BeppoSAX 


2.066 


1.9 ±0.4 


0.7 ±0 


2 


I 


I 


31.4 ±6.8 


310 ± 20 


— 


GRB 010222 


BeppoSAX 


1.477 


1.35 ± 0.06 


1.00 ±0 


06 


I 


I 


94 ± 10 


766 ± 30 





GRB 011121 


BeppoSAX 


0.36 


1.30 ± 0.03 


(1.1) 




II 


II 


9.9 ± 2.2 


793 ± 533 





GRB 011211 


XMM 


2.14 


2.1 ± 0.3 


1.3 ± 


3 


II 


I 


6.3 ± 0.7 


186 ± 24 




CjKB UzU4Uo 


Chandra 


0.69 


1.6 ± 0.9 


1.0 ± 


2 


11 


I 


12.8 ±1.5 
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GRB 030226 


Chandra 


1.98 


2.7 ± 1.6 


0.9 ±0 


3 


II 


I 


14 ± 1.5 


289 ± 66 




GRB 030328 


Chandra 


1.52 


1.6 ±0.3 


1.1 ±0 


2 


II 


I 


43 ±4 


328 ± 55 




GRB 030329 


XMM 


0.168 


(2) 


1.0 ±0 


2 


II 


I 


1.7 ± 0.3 


100 ± 23 





GRB 031203 


XMM 


0.105 


0.5 ± 0.1 


0.8 ±0 


2 


III 





0.010 ±0.004 


158 ± 51 





GRB 050126 


SWIFT 


1.29 


1 1+0.6 


0.7 ±0 


7 


II 





— 


— 


2.34t3;30 


GRB 050223 


SWIFT 


0.5915 


1.4 ±0 


7 


III 




— 


— 


GRB 050315 


SWIFT 


1.949 


"^•'-0.1 


0.91 ±0 


09 


I 


I 


4.9 ± 1.5 


< 89 


4.39 ±0.27 


GRB 050319 


SWIFT 


3.24 




0.96 ±0 


09 


I 


I 


— 


— 


4.67 ±0.27 


GRB 050401 


SWIFT 


2.90 


1 K+O.S 
^•■^-0.2 


1.0 ±0 


3 


I 


I 


41 ± 8 


467 ± 110 


q 0-7+0.24 


GRB 050406 


swift 


2.44 


1 n+'''^ 
J^-^-i.o 


1.0 ± 


4 


II 





— 


— 


2.621?:® 


GRB 050416A 


swift 


0.6535 


0.85 ± 0.05 


1.1 ±0 


1 


II 




0.12 ± 0.02 


25.1 ±4.2 




GRB 050505 


swift 


4.27 


1 t;+0.6 
-'-•^-0.2 


1.0 ±0 


1 


I 


I 






4 Qq+0.48 
4.JM_Q 24 


GRB 050525A 


swift 


0.606 


1.4 ±0.1 


1.1 ±0 


4 


II 


I 


3.39 ±0.17 


127 ± 10 


2.92+?l* 


GRB 050603 


swift 


2.821 


1 a+0.5 


0.7 ±0 


1 


I 


I 


70 ±5 


1333 ± 107 


4.83t?-:?° 


GRB 050730 


swift 


3.968 


r, Y + 0.3 
^- ' -0.2 


0.62 ±0 


05 


I 


I 






4.13 ±0.08 


GRB 050802 


swift 


1.71 


1.7± 0.2 


0.81 ±0 


09 


II 








3.96 ±0.10 


GRB 050814 


swift 


5.3 


0.8 ±0.4 


0.7 ±0 


1 


I 








q Qq+0.56 
'-'■^■^-0.70 


GRB 050824 


swift 


0.83 




0.8 ±0 


2 


II 


II 


0.130 ±0.029 


< 23 


4 00+1.84 
^■''^-0.54 


GRB 050826 


swift 


0.297 


1.13 ±0.04 


1.1±0 


4 


III 










GRB 050908 


swift 


3.344 


1 1+0.3 
J^-J^-0.4 






II 








q q,+0.86 

o.ol_j^ 29 


GRB 050922C 


swift 


2.198 


1.26 ±0.04 


1.3 ±0 


2 


II 


I 


6.1 ± 2.0 


415 ± 111 


2.58 ±0.11 


GRB 051016B 


swift 


0.9364 


n 7+0.2 
^•'-0.3 


0.9 ±0 


2 


II 








q CI +0.64 


GRB 051109A 


swift 


2.346 


1.25 ± 0.07 


1.0 ±0 


2 


I 








o.yo_Q 23 


GRB 051109B 


swift 


0.080 


1.1 ± 0.3 


0.7 ±0 


4 


III 








3 g7+0.27 
-0.40 


GRB 051111 


swift 


1.549 


1.56 ±0.02 


1.1 ±0 


4 


II 


I 








GRB 060108 


swift 


2.03 


1 q+0.3 

^^•■^-0.2 


1.0 ±0 


3 


II 








4.40 ±0.2 


GRB 060115 


swift 


3.53 


1 o+°-2 


1.3 ±0 


3 


I 










GRB 060206 


swift 


4.048 


^■^^-0.06 


0.8 ±0 


2 


I 










GRB 060210 


swift 


3.91 


1-8-0.3 


1.00 ±0 


09 


I 










GRB 060218 


swift 


0.033 


1-^-0.6 


0.5l5±0 


05 


III 








^■'-'^^-0.48 


GRB 060223A 


swift 


4.41 


1.3 ±0.2 


0.9 ±0 


3 


II 








r\ ■7Q~I~0.32 

^- '•^-0.36 


GRB 060418 


swift 


1.489 


1.5 ± 0.2 


0.8 ±0 


9 


II 








3 44+0.14 
•^■^*-0.18 



3. Results 

3.1. Is the clustering still apparent? 

Figure [1] presents the light curves corrected for 
distance effects. As one can clearly see, no clus- 
tering is observed in the early part of the light 
curves. The addition of SWIFT bursts seems 
to have blurred out the clustering properties ob- 
served in paper I and II. If we assume that SWIFT 
is observing GRBs similar to those observed pre- 
viously, the only difference that could explain this 
result is the addition of the early time X-ray 
light curves, since previous observations from Bep- 
poSAX, XMM, Chandra where performed usually 
at least 6h after the GRB alert. 

3.2. Influence of the prompt emission 

SWIFT has shown that the X-ray light curve 
is typically composed by a first steep power- 
law segment, associated with the tail of prompt 
emission, followed by a flat plateau, a second 
steepening (which is the segment usually seen 
before SWIFT) , and a possible late steepen- 
ing related to the j et apert ure (iNousek et al 



20061: lO'Brien et al. 1 ,2 006; Zhang et al. I l2006h . 



Willingale et al. I (l2007f l bave interpreted this be- 
havior with a two-components model. According 
to them, the plateau phase marks the transition 
between a first component, ascribed to the prompt 
emission, and a second one which nature is less 
clear and develops in the afterglow. We have 
retrieved for SWIFT bursts the values of the end- 
time Ta of the plateau phase, and excluded all 
SWIFT data taken before this time. The results 
are displayed in Fig. [2l 

Excluding these data, the two groups reported 
in papers I and II become apparent. We note how- 
ever still some dispersion during the first part of 
the light curves. We interpret this as a conse- 
quence on the error (sometime large) on the Ta 
measurement. Indeed, restricting the sample to 
data collected after Ta + axa (i-e. being very con- 
servative on the Ta value) strongly reduce this dis- 
persion, at the price of a drastic reduction of the 
sample size. In the following, we will use all data 
taken after Ta (i.e. a less conservative hypothe- 
sis) and do not discuss the effects seen in the early 
part of the light curves. We refer to the bright 
group as group I and the dimmer one as group II, 



as in paper II. They are clustering with a mean 
flux of 7.0 X 10-12 erg s'^ cm'^ and 3.1 x lO'^^ 
erg s^i cm^2 for groups / and //respectively, as- 
suming a common redshift of unity (see above) 
and at a time of 1 day. Some bursts do not fol- 
low this relation : GRB 980425, GRB 031203, 
GRB 050223, GRB050826, GRB 051109B, GRB 
060218, GRB 060512, and GRB 060614. They are 
all nearby events, the most distant, GRB 060512, 
being at z = 0.4428 (jBloom et al. II2006I ). Note 



however that GRB 030329 is also at low redshift 
(z — 0.168) but belongs to group II. We consider 
that these low luminosity events form a specific 
group, referred as group III in the following. We 
finally observe two peculiar events. GRB 060605 
belongs to group I but, due to its steep decay, 
changes to group II later on. We note that, con- 
trary to the other bursts, this event is not compat- 
ible with an ISM or wind medium, but is compat- 
ible with a jet effect. This would imply that the 
Boer & Gendre relation is valid before jet effects 
become apparent. GRB 060604 is located between 
group / and group II. On 61 events, this is the only 
one lying within the gap. We note that the light 
curve of this event presents a strong "noise" (i.e. 
small scale temporal variations) during the whole 
observation. A possible explanation to these vari- 
ations would be that this event presents several X- 
ray flares (like e.g. GRB 050904), implying that 
we do not observe the continuum. Under this as- 
sumption, this event would be a group II event 
(and should not be considered for further discus- 
sion). However, we cannot rule out a possible out- 
liers hypothesis and maintain this event within the 
sample. 

The probability that a power law luminosity 
distribution (letting the index be a free parame- 
ter) represents the observed distribution is, at the 
maximum, 8.9 x 10"^^; thus the observed cluster- 
ing in two groups is very significant. To compute 
this probability we impose a lower luminosity limit 
such that the group III bursts are excluded, be- 
cause of the selection effects (see Sec. |4]). 

We checked if this clustering can be related to 
the isotropic burst energy /?iso or the peak en- 
ergy Ep values. Figure [3] presents the Ep distri- 
bution of bursts versus the X-ray afterglow flux 
at one day. One can clearly see that the two 
groups share a similar Ep distribution. This also 
holds for the Eiso distribution (Fig. HJ see also 
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Table 1 — Continued 



Burst name 


X-ray 
satellite 


redshift 


Decay 
index 


Spectral 
index 


X-ray 
group 


Optical 
group 


(10^2 erg) 


Ep 

(koV) 


log(T,) 
(s) 


GRB 060510B 
GRB 060512 


SWIFT 
SWIFT 


4.9 
0.4428 


1.2 ±0.2 


1.7 ±0.4 
0.9 ±0.2 


I 

III 








4 Kc+o.eo 

^•''^-0.54 

3.85+11" 


GRB 060522 


SWIFT 


5.11 


1 o+"-2 

^■^-0.1 


1.1 ±0.2 


I 








2.86l°;« 


GRB 060526 


SWIFT 


3.21 


1 1+0.3 
^■^-0.2 


0.7 ±0.3 


I 








3-»4_Q 33 


GRB 060604 


SWIFT 


2.68 


1.2 ±0.2 


1.3 ±0.2 


7 








4.55 ± 0.29 


GRB 060605 


SWIFT 


3.78 


2 0+°-'^ 
^■'J-o.s 


1.0 ±0.1 


I ? 








4.16 ±0.13 


GRB 060607A 


SWIFT 


3.082 




0.86 ± 0.09 


II 








4.75 ±0.03 


GRB 060614 


SWIFT 


0.125 




0.8 ±0.2 


III 








5.00 ±0.09 


GRB 060707 


SWIFT 


3.43 


0.83 ±0.08 


1.2 ±0.5 


I 








q KO + 0.47 


GRB 060714 
GRB 060729 


SWIFT 
SWIFT 


2.71 
0.54 


1 otr + O.OQ 
-^•^^-0.07 

1.32 ±0.05 


1.4 ±0.3 
1.39 ±0.05 


I 

II 








o 20+0-21 
5.11 ±0.04 




Fig. 2. — Left: the corrected light curves of our sample, using only data taken after the plateau phase for 
SWIFT bursts (see text). The groups reported in papers I and II are now clearly seen. Right: the distribution 
of fluxes observed at one day after the burst. The solid lines are the best Gaussian fit distribution of the 
group fluxes. In the following flgures, red refers to bright events, blue to intermediate luminosity events and 
green to under-luminous events. See electronic version for colors. 
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Sec. ID fo r a discussion on the selectio n effects) 



Recently, Chattopadhvav et al. ( 2007t ) reported 



three classes of GRBs based on their fluence and 
duration. However, their classification and the re- 
sults we report here are not correlated. More- 
over the early SWIFT X-ray light curves associ- 
ated with the tail of the prompt emission do not 
cluster in luminosity. This leads us to consider 
that the afterglow emission is not correlated with 
the prompt emission. 

3.3. Influence of the afterglow properties 

In paper II, we pointed out a possible segrega- 
tion of decay indexes : group I bursts seemed to 
decay faster than group //ones. Thanks to the fast 
monitoring capabilities of SWIFT we can clearly 
rule out this hypothesis. Figure [5] shows that there 
is no difference between the decay index distribu- 
tions for the two groups. Moreover, SWIFT has 
added some events to the group / with a low decay 
(e.g. GRB050315). We find the same results when 
we look at the spectral indexes (Fig. [6|). Clearly, 
the two distributions are similar. In fact, except 
for two of them, all bursts are compatible with 
having a common spectral index of/3 = 1.0±0.2 
(with the flux F oc t-^iy-l^). 



4. Discussion 

4.1. Distance effects and bias 

We present all results using a common distance; 
however, as we are correcting for distance effects, 
one may wonder if this could bias our findings. 

Figure [7] presents the redshift distribution of 
our sample. There is a clear trend : low luminos- 
ity bursts are nearby while bright ones are more 
distant. This is even more obvious if we include 
group ///bursts, which have all a very low redshift. 
This effect can be partly explained by a selection 
effect, as low luminosity events cannot be detected 
at high distance. The lack of dim distant events 
(located in the upper left corner of Fig. [7]) is due 
to that effect. However, we should also have de- 
tected bright nearby events. The nearest group 
I burst, (GRB 991216), is located at z = 1.02. 
While this is alread y lower than the mean S WIFT 
GRB redshift (2.7. [jakobsson et al. Iliooeh . there 
is a clear deficit of bright sources closer z ~ 0.5: 
whether this is related to the nature, or the in- 




S 10" 



10" 



10 



10 10 10 

Time since burst (days) 



Fig. 1. — The light curves of our sample, cor- 
rected for distance effects. We indicate SWIFT, 
XMM-Newton, Chandra and BeppoSAX bursts 
with blue, green, black and red diamonds respec- 
tively. The error bars are not plotted for clarity. 
See electronic version for colors. 
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10"' 10" 
X-ray Flux at one day (10"^' 



10 
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Fig. 3. — Distribution of intrinsic Ep (rest frame) 
values. Blue circles and red squares represent 
group II and group I bursts respectively. The flux 
at one day as been computed by interpolating the 
flux at that date when the burst was observed one 
day after the burst (see text for discussion about 
extrapolation effects), thus lowering the number 
of events available for comparison. 
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Fig. 4. — Distribution of the energy emitted dur- 
ing the prompt event assuming an isotropic emis- 
sion (Eiso)- Symbols are the same as in Fig. [3l 
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Fig. 6. — Distribution of the spectral index of the 
groups / and //. No differences can be seen. Sym- 
bols are the same as in Fig. [31 



3.5r 



3- 



2.5- 



2 



1.5- 



0.5- 




10"' 10" 10' 10' 

X-ray Flux at one day (10"'^ ergs.cm"^. s"') 

Fig. 5. — Distribution of the decay index of the 
groups / and //. No differences can be seen. Sym- 
bols are the same as in Fig. [3l 
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Fig. 7. — Distribution of the redshifts of the 
groups / and //. Symbols are the same as in Fig. 
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trinsic frequency, of group /bursts, is still an open 
question. 

Another important effect to take into account 
is the time dilation and the estimation of the 
distribution parameters. X-ray light curves are 
not taken continuously. Before SWIFT, X-ray 
observations were usu ally performed as snap- 
shots of 30-60 ksec l ong ( De Pasquale et ai~1l2006t 
Gendre et al. I2OO6I) . Because of the time dilation, 



the start and end of each observation were located 
at different times, and the estimation of a param- 
eter at a given time needs extrapolation or inter- 
polation from the value observed at another time. 
However, the decay index can largely vary dur- 



i ng th e observations as shown by lO'Brien et al. 



( 20061 ): this uncertainty can affect the accuracy 
of the distribution parameter estimations. This 
eff ect can be illu s trated by looking at Fig. 6 



of iNardini et al. ( 2006f l , which is equivalent to 



— "---—"I — 

Nardini et al. 


|200d) 


Nardini et al. 


mm 



Fig. [5] of the present work: the results are clearly 
a strong clustering, while 
see only a trend. However 
extrapolate the burst after- 
glow light curves, when no data is available at 
times earlier than 12h in the source rest frame: 
as an example GRB 020405 is observed from 1.00 
day to 1.34 days (source rest frame), GRB 000926 
is observed from 0.69 to 4.36 days (idemj^. The 
net effect of doing extrapolations instead of inter- 
polations is to broaden the distribution, blurring 
in turn the clustering reported here. To avoid this 
effect, we decided, as in paper I and H, to avoid 
extrapolations and to perform only flux interpo- 
lations, when relevant, using the nearest measure- 
ments. This is of course at the expense of the 
sample size, as one can note in Figs. [21 IH El HI [3 
where only half of the total sample is plotted. 

4.2. The optical afterglows 



Kann et al. 1 (I2OO6I) iNardini et cd~l (|2006l) and 



Liang fc Zhang I ( 20061 ) have shown that even op- 
tical afterglow light curves display a clustering ef- 
fect. These authors do not include in their works 
any of our group III bursts: for consistency we 
will restrict the following discussion to groups / 



^ We stress that this example is given only as i llustrative pur- 
pose and do not to questions the results of lNardini et alH 
12006) at optical wavelengths. As the authors noted, opti- 
cal light curves have a better temporal coverage and sam- 
pling, thus in these case they perform only interpolations. 



and //. They found two groups, and we will re- 
fer to them as optical group I (or oT) and opti- 
cal group II (or oil) for the bright and dim group 
respectively. We indicate in Table [1] the optical 
class of each bursts. Several events of our sam- 



ple ar e not incl uded in the work o f INardini et al. 
(|2006l ) or iLiang fc Zhang I (|2006h . and no optical 
classification is r eporte d for them. According to 



iLiang fc Zhang (|2006l) . one day after the burst 
(rest frame), the difference in luminosity between 
groups ol and oil is ~ 26. At the same time, the 
difference in luminosity between the X-ray groups 
is ~ 24. 

From Table m we can see that while an 0/ burst 
can belong to X-ray group / or X-ray group II, oil 
bursts are dim both at optical and X-ray wave- 
lengths. We present in Fig. [5] the distribution of 
the decay index versus the spectral index of these 
three burst classes. Again, no clear separation can 
be observed. 

4.3. The fireball model 

In this section we investigate further on the na- 
ture of the observed cluste rings in the con t ext o f 



the External Shock model. INardini et al. ] (I2OO6I) 
has indicated that around 1 day for almost all 
bursts, we have i/^ < ^optical < Vc < vx', thus 
in the following we will assume this repartition of 
th e characteristic emissio n frequ encies. According 
to Panaitescu fc Kumar I (l2000l) . during the slow 
cooling regime the predicted flux for a fireball ex- 
panding in a uniform interstellar medium (ISM) 



in2.i-i.3pn-2„-(p-i)/2+-3(p-i)/4 

lU 1^28 ^14.6 '-d 



(p+3)/4^1/2 D-1 ^(p+l)/4 



53 



mjy 



(1) 



= 102 



p/2,-(3p-2)/4 



'■^pr)~^u^''' "t 

^IS '^14.6 
(p+2)/4 p-1 ^(p-2)/4 



53 



mJy 



(2) 

while for a fireball expanding in a wind-like 
medium we have: 
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for a fireball expanding in a ISM and 
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2.3-1.2 



■Pn-2„-(P-l)/2 + -(3p-l)/4 
^19, '^14.6 



(p+l)/4 
53 



-^opt, wind 

Fx. wind 



2Ql. 9+0.1)9 

L35 



(3) 



72 V^'' -1/4 -()9-l)/2 p/2 
27 ) d ^opt, 14.6 '^X, 14.6^ 

(10) 



3/4 
B. 53 



for a fireball expanding in a wind-like medium. 



= 1.35 X 102-4-0-8pi5-2^-p/2i-(3p-2)/4 



We extracted from the work of iNardini et al. 



(p+2)/4 p-l ^(p-2)/4 



53 



-1"=B, 



mjy 



(4) 

where D, and t are the distance, observation 
frequency and observation time respectively, and 
p is the electron distribution power law index; we 
assume p = 2.2 for all bursts. We rewrite these 
equations in terms of Eb = cbE and E^ = CgE, 
the energy carried by the magnetic field and rel- 
ativistic electrons respectively {E is the fireball 
total energy): 



iQ2.1+0.7pr) -,, 



-2, -(p-l)/2,-3(p-l)/4 



xE'i+f'E^^;^,El-Pin^/E,,y/^mJy 



Wrn < 



(5) 



F. 



1 n0.4+0.7p n-2, -p/2 ,-(3p-2)/4 
lU 1^28 '^14.6 ''d 



(6) 



1 n2.3+0.8p n-2„-(P-l)/2^-(3p-l)/4 
iU 1^28 '^14.6 ''d 

xE^i^^l"El;l,El~^{A^/E,,)mJy 



(|2006l ) the mean optical fluxes at 4.69 x lO" Hz. 
They are 0.178 mJy and 10.0/xJy for groups ol 
an^ a/7 Respectively. The mean X-ray fluxes (at 
4.8 X 10^^ Hz) are 2.10 ^liy and 0.09 /iJy for groups 
/ and // respectively. 

Combining equation [S] with equation [51 and 
equation [7] with equation [TUl we obtain the re- 
sults listed in Table [2l One can note that the 
constraints on E^. depend weakly on the medium 
density (because of the very low exponent of the 
density parameter). Because i?e < 75 (or Ce < 1), 
we can find a constraint on the total fireball en- 
ergy E (a similar condition is true for Eb^ but this 
is not constraining for E). 



4-3.1. The fireball parameters 

V <Wd);hin each group the fiux variations around 
the mean value are small. This may indicate 
that either Eb^ E^, E, and (or A^, de- 
pending on the surrounding medium) are cor- 
related quantities, or are all fixed to a com- 
mon value. In the former case, because the 
surrounding medium density is fixed before the 
burst occurs, it is this quantity that should be 
considered as the true variable parameter (and 
not the energy injected within the fireball). On 
the other hand, a stellar progenitor produces a 
wind which affects its surrounding mediurn (e.g . 



(^"i < ^HkAiircz -Ruiz et al. 
(7) 



200lUChevalier et al. |[20ol 



Eldridge et al. Il2006ll . Thus, the surrounding den 



1.35 X {^Y^' X 10°-4+0-7Pi?-/^-//2i;(3p-2)/4 



(8) 

We thus find the optical-to- X-ray flux ratio to 



be: 



''opt, ISM 



F- 



X, ISM 



.E 



3/4 



opt, 14.6 "^X, li.G^B, 53 1 p 

\E53 

(9) 



1/2 



sity is (partly) fixed by the progenitor properties. 
As these properties should also affect the energy 
quantity emitted by the fireball (and thus Eb, E^,, 
and E), we find no surprising that Eb, Ee, E, and 
are correlated together. 
If Eb, Ee, E, and (or A,) have fixed values, 
the constraints on the model are strong since this 
restrict to only three types the possible environ- 
ments in which a ORB can occur (according to the 
three observed behaviors). If so, this is a possible 
explanation of why we do not observe a relativistic 
outflow associated with each type Ib/c supernovae 



A^ 
E53, 
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Table 2: Constraints on Eb, E^, E, and the density parameter implied by the clustering observed in X-ray 
and in optical. 



X-ray 
group 


Optical 
group 


Medium 

class 




constraint 




Eb 
constraint 


E 

constraint 


/ 
/ 


ol 
ol 


ISM 
wind 


pi. 2 
pi. 2 


= OME}fnl^'' 

= o.i6eU''aI/'' 


Eb, 53 
Eb, 53 


= 0.88 X 10-\E53/no)^^^ 
= 0.79 X 10-5(£;53/A*)4/3 


0.04 < E53 
0.17 < ^53 


II 
II 


ol 
ol 


ISM 
wind 


pi. 2 

pi. 2 

^e,53 " 


= 0m6E%^^Al^^^ 


Eb, 53 

Eb,53 


= 5.84 X 10-3(^53/no)"/3 
= 4.95 X 10-^{E53/A,)^^^ 


0.002 < £:53 
0.008 < ^53 


II 
II 


oil 
oil 


ISM 
wind 


pi. 2 

pi. 2 

^e,53 - 


= om2Eli\l^^° 

= 0.007 E^i^^Al^^^^ 


Eb, 53 

Eb,53 


= 1.26 X 10-4(£;53/no)2/3 

= 1.07 X io-5(£;53M.)^/3 


0.002 < £^53 
0.009 < E53 



V 2 




0.5 



1 1.5 
spectral index 



Fig. 8. — X-ray decay index versus spectral index. 
In this figure we plot bursts bright both in X-ray 
and optical (red squares), bursts dim both X-ray 
and optical (blue circles), and bursts bright in op- 
tical and dim in X-ray (blue diamonds); see text 
for discussion. 



(jSoderberg et al. ][2003): 

in fact, normal type Ib/c 
supernovae can't develop as a GRB. 

4-3.2. The nature of the clustering 

The flux ratio between the two o ptical groups 

and th e two X-ray groups are similar. P anaitescu fc Kumar 



(j2002[ ) have investigated the surrounding medium 
around 10 GRBs (GRB 970508, GRB 980519, 
GRB 990123, GRB 990510, GRB 991208, GRB 
991216, GRB 000301C, GRB 000418, GRB 000926, 
and GRB 010222). They found that most of them 
can be fit with an ISM (all belong to our opti- 
cal and X-ray bright burst group), while only one 
burst, GRB 970508, requires a wind environment 
(it is an optical and X-ray dim burst). This sug- 
gests that group //-o// bursts are surrounded by a 
stellar wind while group /-o/are located within an 
ISM. Assuming that the only difference between 
these two groups is due to the medium density 
profile (and thus that all other parameters have 
a similar distribution within the two groups), we 
can explain the difference between the two groups 
if we set n^/E^-i ~ 1233AJ (assuming < Vc). 
However in such a case we then cannot explain the 
behavior of group JJ-o / bursts. We also note that 



Panaitescu &: Kumar I (j2002( ) have found that a 
wind medium can also accurately describe the sur- 
rounding environment of GRB 991216 and GRB 
010222 (two group I-oI bursts). Thus, the sur- 
rounding medium type cannot explain the origin 
of the observed clustering. 

In X-ray, the flux is mostly dependent on Ee 
and E (see eqn. [5] and [5]). If one assumes that Eiso 
is a good estimator of E, then since the observed 
distributions of Eiso are similar for the two groups 
(see Fig. 0]), different values of E cannot explain 
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the differences. In turn this impHes that varies 
within the two groups (see Table [2|). This however 
again cannot explain the origin of group II- ol: if 
the clustering is due only to a difference in the 
values, then we cannot expect from eqn. [5] and [HI 
and [7] and [H that a dim X-ray afterglow can be 
bright in optical. The constraints listed in Table 
[2] indicate that the difference between groups I-oI 
and II-oII on one hand, and group II-oI on the 
other hand can be also due to a different value of 
Eb- 

We thus propose that the presence of different 
groups is ascribed to different families of E^ and 
Eb values: 

• a family of 'magnetized fireball' that pro- 
duces the group II- ol. In such a case the 
fireball transfers only a low fraction of its 
energy into relativistic electrons. 

• a family of 'less magnetized fireball' that 
produces the groups I-oI and II-oII. The 
fraction of total energy going into magnetic 
fields is roughly one order of magnitude 
lower than group II- ol. These two groups 
can be related to an high and low fraction 
of energy going in relativistic electrons, re- 
spectively. 

Interestingly, the group II-oII is not very nu- 
merous (compared to the other two groups), so 
that most of GRBs of our sample are then either 
'magnetized but not electron-energized' or 'not 
magnetized but electron-energized'. 

5. Applications 

5.1. The estimation of GRB source red- 
shifts 

The clustering observed in X-ray can be used 
for distance estimation. To date, most of the 
redshift measurements made on GRB afterglows 
were done by optical spectroscopic or photometric 
observations. However, because of the Lymann 
alpha cut-off, this method cannot be applied to 
high redshift events {z >^ 5 — 7), unless doing 
IR spectroscopy. Moreover, not all GRBs can be 
followed in optical (e.g. t he so-called dark-bursts. 



De Pasauale et al. 2003) , either because they are 



too faint at the time of the observation (or even 



dark), or because there is no large enough tele- 
scope available, or because of the position of the 
event respective to the Sun. Hence many GRBs 
have no known distance, stressing the need for in- 
trinsic distance indicators. After the launch of the 
SWIFT satellite, thanks to its fast re-pointing ca- 
pabilities, nearly all GRBs have an homogeneous 
X-ray follow-up. Hence, a redshift measurement 
method based solely on X-ray observations could 
be very interesting if one wants to use a large 
sample of GRB sources for cosmological studies. 
We propose here to use the reported clustering of 
GRBs X-ray light curves in three classes to build 
a redshift estimator. 

The method is based on the redshift needed 
for a burst to belong to one (or both) of the two 
groups. The steps are easy to follow: we compute 
from the X-ray observations (i.e. obtained from 
SWIFT) the flux at 1 day (observer frame) in the 
2-10 keV band, either using a mean spectral index 
of ^ 1.2, or the exact one obtained from X-ray 
spectral fitting. We note that even if the spec- 
tral index is not known, this has little infiuence on 
the actual fiux at one day (within reasonable lim- 
its). Table [31 which provides the redshift needed 
to comply with the relation for both groups for a 
given observed flux, is directly usable to estimate 
the redshift. 

We calibrated this method by deriving the es- 
timated redshift for the bursts of our sample (for 
which the group is known), and comparing this 
value with the measured redshift. The results are 
displayed in Fig. [9l The estimated redshift agrees 
with the measured one for most of the bursts. The 
only discrepancies arises at low redshift. In Sect. 
U] we already noted that group III bursts are lo- 
cated at small distances (z<0.5) and do not follow 
the same relation as groups I and II: as a con- 
servative approach we prefer to restrict the valid- 
ity of our method to source located at redshifts 
larger than 0.5. We list in Table [3] the observed 
fiux one day after the burst (observer frame) and 
the associated redshift needed for a burst to com- 
ply with the observed clustering for group I and 
II sources. Because of the limitation to redshifts 
larger than 0.5, all bursts brighter than ^ 2x 10^^^ 
erg s~^ cm~^ belong to group I. However, the be- 
longing group is not correlated with easily observ- 
able quantities (absorption, peak energy, isotropic 
energy, spectral or decay indexes), so it can be 
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fixed only through broad-band modeHng. In fact 
this leads to two redshift estimates. If the bursts 
has an optical afterglow there is a possible way to 
decide which of the estimate is valid: a burst with 
flux of ~ 2 X 10~^^ erg s~^ cm~^ which belongs 
to group / has a redshift of ^ 2. At that distance, 
the Lymann a cut-off appears to be at 360nm, in 
the UVOT observation bands. Hence, if a GRB 
afterglow is observed by UVOT (or in the B band 
for ground based telescopes) then it cannot be- 
long to group I, and the ambiguity on the redshift 
determination is cleared. 

5.2. Estimated redshift of pre-SWIFT 
bursts with unknown distance 



We retrieved from Gendre et al. I (l2006l) and 



De Pasquale et al~1 ( 20061 ) the light curves of sev- 



eral bursts with good temporal sampling and spec- 
tral informations (table |3|), and estimated their 
redshift assuming they belong either to group I or 
to group II. For only one burst we compute an 
associated redshift lower than 0.5 (GRB 020410), 
and thus we do not consider it in the following 
discussion. 

The estimated redshift of GRB 980519 is either 
3.8ib 0.7 or 1.4 zfc 0.2. This burst was observed in U 
band [jaunsen et aL 



( 200lh : lKann et al. 



(l2006l ). 

As the Lymann a cut-off cross the U band at 
z ^ 2.8, this burst should not have been observed 
in the high distance hypothesis. Thus, U observa- 
tions can indeed solve the problem of group clas- 
sification, and we propose a redshift measurement 
for GRB 980519 of 1.4 ± 0.2 based on X-ray obser- 
vations. For the same reason we argue that GRB 
040827 is at = 1.9 ± 0.3 and not z = 8 ± 2, since 
an optical afterglow h as been observed in optical 



( Malesani et al 1120041 ) 



GRB 001025A was observed by XMM-Newton 
but no optical afterglow was dete cted. This event 
was t hus classified as a dark burst (jPedersen et al. 
20061 ). We note that this one has a large red- 
shift value, most of all if it belong to group 
I. In such a case, with a calculated redshift 
of 5.8 ± 0.8, we can explain the classification 
as dark burst by the Lyman alpha cut-off at 
^ 6200A. However, since we cannot draw con- 
clusions from the absence of detection at optical 
wavelengths, and as dark bursts can also b e in- 

2003h . 



Table 3: Flux to redshift conversion. This table 
should be used as an estimate of the redshift for 
bursts observed by any X-ray observatory. The 
flux is given in the 2-10 keV band at 1 day after 
the burst (observer frame). The redshift has been 
calculated for an energy index of 1.2, the uncer- 
tainty is 30%. 



Flux 


Group I redshift 


group II redshift 


(erg s^^ cm^^) 






1 X 10-^^ 


— 


4.43 


2 X 10-1'' 


— 


3.28 


3 X 10-1'' 


— 


2.72 


4 X 10-1'' 


— 


2.35 


5 X 10-1" 


— 


2.12 


6 X lO-i-i 


7.80 


1.96 


7 X 10-1" 


7.09 


1.83 


8 X 10-1" 


6.38 


1.75 


9 X IQ-i" 


6.05 


1.68 


1 X 10-13 


5.78 


1.62 


2 X 10-13 


4.10 


1.26 


3 X 10-13 


3.48 


1.09 


4 X 10-13 


3.05 


0.96 


5 X 10-13 


2.82 


0.89 


6 X 10-13 


2.59 


0.82 


7 X 10-13 


2.42 


0.77 


8 X 10-13 


2.29 


0.73 


9 X 10-13 


2.10 


0.69 


1 X 10-12 


1.99 


0.66 


2 X 10-12 


1.50 


0.50 


3 X 10-12 


1.29 




4 X 10-12 


1.16 




5 X 10-12 


1.06 




6 X 10-12 


0.99 




7 X 10-12 


0.94 




8 X 10-12 


0.89 




9 X 10-12 


0.83 




1 X 10-11 


0.79 




2 X 10-11 


0.56 




3 X 10-11 


0.50 





trinsically fainter (jPe Pasquale et al 



we 



cannot exclude the hypothesis that this event is 
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Fig. 9. — Comparison of the estimated redshift 
versus the measured one for the bursts defining 
our sample. The solid line corresponds to equality 



Table 4: Redshift estimates derived from the rela- 
tion using either the group I and II hypotheses for 
pre-SWIFT bursts without known distance. 



GRB name 


Redshift 


estimate 




group I 


group II 



GRB 980329 


4.2 ± 1.2 


1 


2 ± 





2 


GRB 980519 


3.8 ± 0.7 


1 


4 ± 





2 


GRB 990704 


3.5 ± 0.9 


1 


3 ± 





3 


GRB 990806 


4 7+1-6 


1 


6 ± 





3 


GRB 001109 


2.3 ± 0.7 





8 ± 





2 


GRB 001025A 


5.8 ± 1.8 


2 


2 ± 





4 


GRB 020322 


5.0 ± 1.5 


1 


5 ± 





3 


GRB 020410 


0.5 ± 0.4 




< 


1 




GRB 040106 


3.4 ± 0.5 


1 


± 





2 


GRB 040223 


5 5+2.0 


1 


7 ± 





2 


GRB 040827 


8.0 ± 2.0 


1 


9 ± 





3 



located at z = 2.2 ± 0.4, hence belonging to group 
11. 

Finally, the mean redshift for SWIFT bursts is 
2.3, while the pre-S WIFT mean redshift was 1.2 
()GruDe et al. II2007I) . Lets assume that all bursts 
listed in Table [4] are group //bursts, a very con- 
servative hypothesis, the mean redshift, in this 
case, is 1.5 : this is already larger than the pre- 
SWIFT mean redshift. Since it is reasonable to 
consider that at least some of these events belong 
to group /, this should be considered as a lower 
limit. Before SWIFT, the X-ray/optical follow-up 
was late, and only the brightest events were ob- 
served. As more distant events will appear fainter, 
this introduced a bias against distant events. In- 
deed, the fact that these bursts without known 
redshift are slightly more distant indicates that 
the pre-SWIFT redshift distribution was biased 
toward low redshift due to selection effects, and 
that SWIFT cleared this bias thanks to its fast 
follow-up. 

6. Conclusions 

We have investigated the clustering of afterglow 
light curves observed previously (papers I and II) 
in X-ray and in optical from BeppoSAX, XMM- 
Newton and Chandra data, using the newly avail- 
able, well-monitored SWIFT light curves. Adding 
SWIFT bursts to the previous sample reported in 
paper II, we still confirm our previous findings. 
On a sample of 61 events the X-ray light curves 
cluster in two groups, with a sig nificance larger 
than 6 a. IWillingale et al. I (|2007n describe the X- 
ray light curve as a sum of two components, the 
first one being related to the prompt emission: we 
found that it is the second component that clus- 
ters in luminosity. This finding is also supported 
by the fact that the clustering is not related to the 
properties of the prompt emission. 



Table 5: Redshift estimates derived from the rela- 
tion using either the group I and II hypotheses for 
SWIFT bursts with known distance. 



GRB name 


Group I 


Group II 


Measured 




estimate 


estimate 


redshift 


GRB 070529 


7.8 


1.96 


2.44 


GRB 070611 


6.05 


1.68 


2.04 


GRB 070721B 




3.28 


3.62 
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A similar clustering was observed also at op- 
tical wavelengths. We compared the classifica- 
tion within each group in X-ray and in optical, 
and found three classes: bright optical and X-ray 
afterglows, dim ones, and optically bright-H.-Tay 
dim ones. We propose that this clustering is re- 
lated to the fireball properties and the surround- 
ing medium density, which can be either constant 
or correlated. We used the mean X-ray and op- 
tical fluxes of each group to put constraints on 
the values of es and eg, and expressed them in 
terms of the total fireball energy and the surround- 
ing medium density. We stress that the exten- 
sion of this work at low frequency (radio, sub- 
millimeter and infrared) may help solving the ex- 
act origin of the clustering by strongly constrain- 
ing the medium density and the position of the 
synchrotron self absorption and injection frequen- 
cies. 

Using the observed X-ray afterglow properties 
of GRBs, we propose a new, simple, method for 
the determination of the source redshift based on 
X-ray data; optical photometry in U and V bands 
may help to clear the degeneracy between the two 
estimates found. This method has been estab- 
lished on a sample of bursts with known redshift, 
and we find an excellent correlation. We apply 
this method to a sample of GRB source of un- 
known redshift. We propose an estimation of the 
redshift for GRB 980519 (1.4 ± 0.2) and for GRB 
040827 (1.9 ±0.3). 

While we were building this distance indicator 
and writing this paper, several GRB were detected 
by SWIFT and their distance derived using opti- 
cal spectroscopy. We thus have a sample of bursts 
which allows to test this correlation independently. 
To simulate actual conditions for an unknown 
event we used the SWIFT c ount rate light curv e 



available from their web site (jEvans et al. II2007I ). 



and we apply a mean spectral index of 1.2 for 
the count-to-flux conversion. The results are indi- 
cated in Table [5l indicating an agreement between 
10.3 and 24.5 % (in accordance with the error esti- 
mated in 15. 21 we estimate far superior than other 
methods e.g. lAtteia ii2003i) . 
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